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We deal with randomness-quantifiers and concentrate on their abihty do discern the hallmark of 
chaos in time-series used in connection with pseudo random number generators (PRNG). Workers 
in the field are motivated to use chaotic maps for generating PRNGs because of the simplicity of 
their implementation. Although there exist very efficient general-purpose benchmarks for testing 
PRNGs, we feel that the analysis provided here sheds additional didactic light on the importance 
of the main statistical characteristics of a chaotic map, namely, i) its invariant measure and ii) 
the mixing constant. This is of help in answering two questions that arise in applications, that 
is, (1) which is the best PRNG among the available ones? and (2) If a given PRNG turns out 
not to be good enough and a randomization procedure must still be applied to it, which is the 
best applicable randomization procedure?. Our answer provides a comparative analysis of several 
quantifiers advanced in the extant literature. 



I. INTRODUCTION 

Chaos theory started more than thirty years ago and 
changed our world view regarding the role of random- 
ness and determinism. As the statistical characteristics 
of chaotic systems were better understood [l, |2„ ^ a 
wide variety of situations emerged in which chaos, instead 
of stochastic systems, became a "controller of noise" . 

Chaos illustrates the rather striking fact that com- 
plex behavior arises from simple rules when nonlinear- 
ities are present. Since simple chaotic maps are capa- 
ble to generate stochastic-like signals, implementations 
based on chaotic systems are usually less involved than 
those based in more complex algorithms [sl, @, 0| • One 
tryes to apply this notion to generate PRNGs because 
random numbers are widely used not only in cryptog- 
raphy and Monte Carlo applications but in less obvious 
applications [I,[l,[l3i[iiijl2]- We mention just a couple of 
them: 1) in spread spectrum techniques, a binary signal 
is mixed with a random number sequence to spread the 
spectrum over a wider frequency range. Using different 
random number sequences it is possible to share a com- 
munication channel among several users [l^ Hj [S ■ 
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Reduction of electromagnetic interference is another im- 
portant benefit of the spread spectrum effect [13, ; 2) 
Consider a low frequency signal immersed in a high fre- 
quency digital noise. Sampling at time intervals defined 
by a random number sequence, the resultant signal be- 
comes filtered without using any coil or capacitor that 
are expensive, specially in power systems |l9{ . 

Truly random numbers are not attainable from com- 
puters and it is likely that we will ever be able to get 
them from "natural" sources, since one commonly as- 
sumes that any system is governed by underlying phys- 
ical rules and consequently it is deterministic. A suc- 
cessful strategy to build up a PRNG is to start with 
the time series of a simple nonlinear chaotic map and 
to apply to it an adequate randomizing procedure so as 
to "heighten/boost" its stochastic nature. Such strat- 
egy requires a quantitative evaluation of the improvement 
achieved after effecting the procedure. In [iO] the Statis- 
tical Complexity Measure originally proposed by Lopez 
Ruiz et al. [2l| and later modified by Lamberti et al. [l^l 
was used to quantify the effectiveness of such randomiz- 
ing modus operandi when applied to a Lorenzian chaotic 
system. It was also shown there that a widely employed 
course of action -the mixing of two chaotic signals- is 
not effective in this respect, contrary to what one might 
expect. In this vein it is important to note that general- 
purpose tests available in the open literature (23j are not 
designed taking into account the particular characteris- 
tics of a chaotic map. Instead, one can appreciate in [l^l 
the fact that the deterministic nature of chaotic dynam- 
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ics leaves special manifestations in the associated time 
series that can be revealed only by recourse to adequate 
statistical quantifiers. 

In [255 chaotic maps were randomized by means of two 
different randomizing procedures, namely, discretization 
and skipping. The idea of concocting an "information" 
plane, called the entropy-complexity plane, was advanced 
in order to use it as a means to ascertain the effectiveness 
of each of these two modus operandis. The main differ- 
ence of this information plane with other Complexity- 
Entropy diagrams is the joint use in it of two different 
probability distribution functions, both associated to the 
pertinent time series. 

Other important tools at our disposal are to be men- 
tioned as well. In a recent and excellent report, Marwan 
et al. reviewed applications of so-called Recurrence Plots 
for a wide variety of fields and endeavors. They also pro- 
posed several measures to quantify the recurrence plots' 
characteristics 26] . Additionally, two useful information- 
theoretic quantifiers of randomness, the rate entropy and 
the Excess Entropy were proposed in [27| as coordinates 
of a Complexity-Entropy diagram. 

In the present work we explore combinations of all the 
above mentioned quantifiers with the purpose to answer 
the following questions: 1) among several chaotic maps, 
just which is the one that generates the best time se- 
ries?; 2) which is the best strategy -Discretization or 
Skipping- to randomize a given chaotic time series?. 
The ensuing quantifiers-testing will be made by means 
of two representative chaotic maps (and their iterates). 

The paper is organized as follows: the statistical prop- 
erties of a chaotic map and its iterates is reviewed in 
section |lll Section IIIII describes each of the analyzed 
quantifiers. Section llVI deals with results for two repre- 
sentative maps and, finally, conclusions are presented in 
Section IVl 



II. STATISTICAL PROPERTIES OF A 
CHAOTIC MAP 

Let / be a chaotic map on the interval [0, 1]. Suppose 
the map has an invariant measure ^i{x). Then the map 
is ergodic if for any integrable test function Q{x), and 
for an arbitrary initial condition xq (up to a set of zero 
/.t-measure) , the time average is equal to the ensemble 
average: 



(Q) 



(1) 



Equation ([T]) is a consequence of the famous Birkhoff er- 
godic theorem (28j . Mixing is an even stronger require- 
ment than ergodicity. A map is called "mixing" if any 
smooth initial probability density p{x) converges to the 
invariant measure ii{x) after enough successive iterations. 
Mixing implies ergodicity. The reverse, however, is not 
true ■ 

There exists an equivalent definition of mixing via cor- 
relation functions. Let (t>i{x) and 4>2{x) be two integrable 



test-functions and define the generalized correlation func- 
tion of the map / by 
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where 
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The map is mixing if, for arbitrary 0i and (^2, 



lim C((?!)i, 02, = 



(2) 
(3) 

(4) 



Let us stress that it is not easy to prove that / is a mix- 
ing map because the mixing condition given in Eq. (j4|) 
must be fulfilled for arbitrary test functions. Formally, 
every mixing map / has an associated Perron- Frobenius 
operator C [29| that determines the time evolution of any 
initial density Pq{x) towards the invariant measure /x(a;): 



Pn+l = ^[Pn] 



(5) 



The explicit formal expression for the Perron-Frobenius 
operator for a one-dimensional map / is given by [1^ 



E 



[Po{x)] 

\nx)\ 



(6) 



This operator C has a set of eigenf unctions ^a{x) and 
eigenvalues r^a- The invariant measure p{x) is the eigen- 
function corresponding to the largest eigenvalue 770 = 1. 
The full set of eigenfunctions and eigenvalues may be 
used as a basis to express any density: 

Po(a;) = ^Ca^a{x) , 

Pn{x) = £"po(a;) = ^?72cq = coV'o(a^) + -Rn • (7) 

a 

The eigenvalue with the second largest absolute value, 771 , 
has a " distinguished" physical meaning: it is related with 
the mixing constant rmix that governs the relaxation of 
exponentially mixing maps: 



l^ril ^ hil" ~ exp(-n/r„„) . 



(8) 



The invariant measure p.{x) gives the histogram of the 
time series and the ideal PRNG must have p.{x) = const. 
The mixing constant rmix gives the transient characteris- 
tic time Il4l . ,[11, [33, Ell and its ideal value is rmix = 0. 
In many applications of PRNGs both the invariant mea- 
sure and the mixing constant are relevant. 

The analytical expression of the invariant measure p{x) 
of a given map / is usually not known. Exceptions are 
the logistic map in full chaos, and the piecewise-linear 
maps. The mixing constant r^nix has been analytically 
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obtained only for piecewise-linear maps. For other maps 
it must be numerically obtained by means of a piecewise 
linear approximation of the map [ll, 0| ■ 

It is then obviously convenient to have quantifiers for 
measuring the uniformity of the invariant measure /i(a;), 
and the mixing constant r„iix, for any chaotic map. 
These quantifiers are useful to compare time series com- 
ing from different chaotic maps and also to assess the 
improvements produced by randomization procedures. 

It is possible to show that the invariant measure of 
is identical to the invariant measure of /. Also, the mix- 
ing constant rmix for J'' is lower that the mixing constant 
for /. The iteration of a map is one of the randomiza- 
tion procedures proposed in the literature, being used to 



diminish 



This procedure is also known as Skip- 



ping because iterating is tantamount to skipping values 
in the original time series, which does not change fJ-{x) 
and, consequently, is not conductive to a randomization 
of chaotic maps with fj.{x) ^ const. In this paper we will 
use "skipping" as a method of quantifier-analysis. 



III. QUANTIFIERS FOR THE INVARIANT 
MEASURE AND MIXING CONSTANT 

In this section we review several quantifiers proposed 
for measuring the main statistical properties of chaotic 
PRNGs. The quantifiers are classified according to their 
origin into three classes: (1) quantifiers based on Infor- 
mation Theory f2ll, [2^, ; (2) quantifiers based on Re- 
currence Plots [26i, i32i] ; (3) quantifiers based on intrinsic 
computation [27| . 



A. Quantifiers based on Information Theory 

They are appropriate functionals of the probability dis- 
tribution function (PDF). Let {xi} be the time series 
under analysis, with length M . There are infinite possi- 
bilities to assign a PDF to a given time series, a subject 
that will be given due consideration below. In the mean- 
time, suppose that the PDF is discrete and is given by 
P = {pi'ji = • • ■ ,N}. One defines then various quanti- 
ties, namely, 

1. Normalized Shannon Entropy H[P]. Let S[P] be 
the Shannon Entropy 



N 



s[p] = -J2p^ Hp^ 



(9) 



Is is well known that the maximum Smax = In(A^) 
is obtained for P^ — {l/N, ■ ■ ■ , that is, the 

uniform PDF. A "normalized" entropy H[P] can 
also be defined in the fashion 



2. Statistical Complexity Measure. A full discus- 
sion about Statistical Complexity Measures exceeds 
the scope of this presentation. For a comparison 
amongst different complexity measures see the ex- 
cellent paper by Wackerbauer et al. . In this pa- 
per we adopt the definition by Lopez Ruiz-Mancini- 
Calbet seminal paper [2l| with the modifications 
advanced in so as to ensure that the concomi- 
tant SCM- version becomes (i) able to grasp essen- 
tial details of the dynamics, (ii) an intensive quan- 
tity and, (Hi) capable of discerning both among 
different degrees of periodicity and chaos [g^l . The 
ensuing measure, to be referred to as the inten- 
sive statistical complexity, is a functional C[P] that 
reads 



C[P] = Qj[P,P,]-H[P], 



(11) 



where Q j is the "disequilibrium" , defined in terms 
of the so-called extensive Jensen-Shannon diver- 
gence (which induces a squared metric) ^2^. One 
has 

Qj[P,P,] ^ Qo-{S[{P+Pe)/2]-S[P]/2-S[Pe]/2}, (12) 

with Qo a normalization constant (0 < Qj < 1) 
that reads 



= -2 



N + 1 
N 



ln(iV + 1) - 21n(27V) +liiN 



H[P] = S[P]/S„ 



(10) 



(13) 

We see that the disequilibrium Qj is an intensive 
quantity that reflects on the systems's "architec- 
ture" , being different from zero only if there exist 
"privileged" , or "more likely" states among the ac- 
cessible ones. C[P] quantifies the presence of cor- 
relational structures as well [13, [SJ • The opposite 
extremes of perfect order and maximal randomness 
possess no structure to speak of and, as a conse- 
quence, C[P] = 0. In between these two special 
instances a wide range of possible degrees of physi- 
cal structure exist, degrees that should be reflected 
in the features of the underlying probability distri- 
bution. In the case of a PRNG the "ideal" values 
are H[P] = 1 and C[P] = 0. 

As pointed out above, P itself is not a uniquely defined 
object and several approaches have been employed in the 
literature so as to "extract" P from the given time series. 
Just to mention some frequently used extraction proce- 
dures: a) time series histogram |35|, b ) binary symbolic- 
dynamics l36[ , c) Fourier analysis 33, rfj wavelet trans- 
form [ssl. I39II , e) partition entropies|40| . f) permutation 
entropy [4ll |42|. g) discrete entropies [ij], etc. There is 
ample liberty to choose among them. In [25| two prob- 
ability distribution were proposed as relevant for testing 
the uniformity of fJ-{x) and the mixing constant: (a) a 
P based on time series' histograms and (b) a P based 
on ordinal patterns (permutation ordering) that derives 
from using the Bandt-Pompe method [41] . 
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For extracting P via the histogram divide the interval 
[0, 1] into a finite number A'^^m of non overlapping subin- 
tervals A,: [0,1] = U,=T and A^f]Aj = Vi 7^ j. 
Note that N in eq. Q is equal to Nbm- Of course, in 
this approach the temporal order of the time-series plays 
no role at all. The quantifiers obtained via the ensuing 
PDF are called in this paper i/C*'"*) and CC**"*'. Let 
us stress that for time series within a finite alphabet it 
is relevant to consider an optimal value of Ni,in (see i.e. 

In extracting P by recourse to the Bandt-Pompe 
method the resulting probability distribution P is based 
on the details of the attractor-reconstruction proce- 
dure. Causal information is, consequently, duly in- 
corporated into the construction-process that yields P. 
The quantifiers obtained via the ensuing PDF are called 
in this paper H'-^P^ and C'^^P\ A notable Bandt- 
Pompe result consists in getting a clear improvement 
in the quality of information theory-based quantifiers 

[H H S, Si, 113, S, m, [131 . 

The extracting procedure is as follows. For the time- 
series {xt : t = 1, ■ ■ ■ , M} and an embedding dimension 
D > 1, one looks for "ordinal patterns" of order D [4ll, 
IHj [m generated by 



(s) 1-^ {xs_i^D_i), 



s-(D-2) 



(14) 



which assign to each "time s" a ZJ-dimensional vector of 
values pertaining to the times s,s — l,---,s — (£* — 1). 
Clearly, the greater the Z?— value, the more information 
on "the past" is incorporated into these vectors. By the 
"ordinal pattern" related to the time (s) we mean the 
permutation tt = (ro, ri, ■ • ■ , ro-i) of (0, 1, • ■ • , D — 1) 
defined by 



(15) 



In order to get a unique result we consider that < 
ri_i if Xs-n = Xs^n-i- Thus, for all the D\ possible 
permutations tt of order D, the probability distribution 
P — {p{tt)} is defined by 

^{s\s<M-D + l;{s) has type ^} 
= M-D + l ■ 

In the last expression the symbol ft stands for "number" . 

The advantages of the Bandt-Pompe method reside 
in a) its simplicity, b) the associated extremely fast 
calculation-process, c) its robustness in presence of ob- 
servational and dynamical noise, and d) its invariance 
with respect to nonlinear monotonous transformations. 
The Bandt-Pompes methodology is not restricted to time 
series representative of low dimensional dynamical sys- 
tems but can be applied to any type of time series (reg- 
ular, chaotic, noisy, or reality based), with a very weak 
stationary assumption (for k = E> the probability for 
Xt < Xt+k should not depend on t |4l|). One also assumes 
that enough data are available for a correct attractor- 
reconstruction. Of course, the embedding dimension D 



plays an important role in the evaluation of the appro- 
priate probability distribution because D determines the 
number of accessible states Dl. Also, it conditions the 
minimum acceptable length AI ^ D\ of the time series 
that one needs in order to work with a reliable statistics. 
In relation to this last point Bandt and Pompe suggest, 
for practical purposes, to work with 3 < D < 7 with a 
time lag t = 1. This is what we do here (in the present 
work D = 6 is used). 



B. Quantifiers based on Recurrence Plots: 

Recurrence Plots were introduced by Eckmann et al. 
[S^ l so as to to visualize the recurrence of states dur- 
ing phase space-evolution. The recurrence plot is a two 
dimensional representation in which both axes are time- 
ones. The recurrence of a state appearing at two given 
times ti, tj is pictured in the two-dimensional graph by 
means of either black or white dots, where a black dot 
signals a recurrence. Of course only periodic continuous 
systems will have exact recurrences. In any other case 
one detects only approximate recurrences, up to an error 
e. The so-called recurrence function can be mathemati- 
cally expressed as: 



R(*,j)-e(£-||^(z)-i?(j-)||), 



(17) 



with ~x{i) e 5R™ and i,j — 1, - ■ ■ ,N. Being N the num- 
ber of discrete states ~x{i) considered, || • || is a norm, 
and 8(») is the Heaviside step function. 

In the particular case of the PRNGs analyzed in this 
paper only IV series are considered but the recurrence 
function-idea can be extended to I?-dimensional phase 
spaces or even to suitably reconstructed embedding phase 
spaces. Of course, the visual impact produced by the 
recurrence plot is insufficient to compare the quality of 
different PRNGs, because of the "small scale" structures 
that may be present in our scenario. Several kind of 
measures have been defined to quantify these small scale 
structures [26], each measure being a functional of the 
recurrence function (Eq. (jl7p ). In this paper two kind of 
Recurrence Plot-measures are considered, namely, 

1. Measures based on the recurrence density (mea- 
sured by the number of points in the Recurrence 
Plot). In this paper we use the Recurrence Rate 
{RR), given by: 



1 ^ 



(18) 



Note that in the limit — » 00, RR is the proba- 
bility that a state recurs to its £-neighborhood in 
phase space. For PRNGs the ideal value would 
be RR = 0. But in practice, if no points are to 
be found in the Recurrence Plot, a larger e must 
be adopted in order that the quantifier may make 
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2. Diagonal Measures: these are measures related to 
the histogram P{e, I) of diagonal line lengths, given 
by 



N 



Pie, I) = ^ [1 - R^-l,J-l{e)] [1 - R^+l,M£)] 



/-I 



l[R,+k^,+k{e) . (19) 



A:=0 

Processes with uncorrelated or weakly correlated 
behavior originate no (or just very short) diagonals, 
whereas deterministic processes give rise to "long" 
diagonals and smaller amount of single, isolated re- 
currence points. In this paper three measures based 
on the statistics of diagonal lines are considered: 

(a) the deterministic quantifier DET, the ratio 
of recurrence points that form diagonal struc- 
tures of at least length to all recurrence 
points 



DET = 



(b) The average diagonal line length L 
El: Pie J) 



L = 



(c) The entropy ENTR given by 



AT 



ENTR =- P(e, /) In P(e, /) 



(20) 



(21) 



(22) 



l = lr, 



C. Quantifiers based on intrinsic computation 

We consider in this paper two quantifiers introduced in 
[131, i.e., the entropy rate and the entropy excess E. 
They are defined for time series with a finite alphabet A, 
which is not a limitation because the cc^'s may be thought 
of as real numbers only in analytical studies. In any 
practical case they are in fact floating point numbers and, 
consequently, they have only a finite number of allowed 
v4— values. A subsequence — {xi, Xi+i, • • • ^xi+l} is 
called an L-block. Let P{s^) denotes the probability of 
a particular L-block. Then the block entropy H{L) is: 



(23) 



The sum runs over all possible blocks of length L > 
and H{0) = by definition. For stationary processes 
and sufficiently large L, H{L) ^ L. On the other hand, 
the entropy rate h^^ is defined as 



H{L) 

L ' 
lim ft,^ (L) 



(24) 



The entropy rate is also known as the metric entropy 
in dynamical systems' theory and it is equivalent to the 
thermodynamic entropy density familiar from equilibrium 
statistical mechanics 27] . The entropy rate provides a re- 
liable and well understood measure of the randomness or 
disorder intrinsic to a process. However, this tell us lit- 
tle about the underlying system's organization, structure 
or correlations. A measure of system's organization may 
be obtained by looking at the manner in which h^{L) 
converges to its asymptotic value h^. When only obser- 
vations over length L-blocks are taken into account, a 
process appears to have an entropy rate of h^{L) larger 
than the asymptotic value of /i^. As a result, the process 
seems to be of a more random nature than it really is by 
the "excess" amount of h^{L) — bits. Summing these 
entropy over-estimates over L one obtains the excess en- 
tropy [52| : 



E 



(25) 



L=l 



D. Expected behavior for PRNG 

Summing up, the quantifiers to be compared here 
are: ijC'*^*), c'^hist)^ fjiBP)^ C'^^p) ^ jj^j^^ ENTR, 

L, hfj^, and E. These quantities should tell us how 
good our PRNG is as compared to the ideal condition 
n{x) = const, rmix ~ 0. jji^ist) ^j^g natural quan- 
tifier to measure a non constant ^{x), with value 1 for 
the ideal PRNG. It does not depend on the order of ap- 
pearance of a given time-series event, but only in the 
number of times such event takes place. As for if ^ 
it is not able to uncover any change in rmia; -values. Thus, 
to get a good representation plane we ought to demand a 
quantifier that changes with rmix and not with pl(x). To 
look for such a kind of quantifier we must study i/C""*)^ 
(j{Mst)^ jj{BP)^ C(BP)^ JIH^ DET, ENTR, L, h^, and E 
as functions of rmix ■ A family of iterated maps may 
be used to that end because they share the same invari- 
ant measure and rmix is a decreasing function of d. The 
best quantifier for r,nix would be that which has maximal 
variation over the entire family of maps. 



IV. APPLICATION TO LOGISTIC MAP AND 
THREE WAY BERNOULLI MAP 

In this Section we present results for the families of it- 
erates of two chaotic well-known maps, the Logistic Map 
(LOG) and the Three Way BernouUi Map (TWB), that 
have been selected, among other possibilities, because 
they are representative of two different classes of systems. 



• LOG is given by: 



''Tl + l 



4 Xyi (1 Xyi) , 



(26) 
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d TWB LOG 

1 0.56789 0.333333333 

2 0.31848 0.111111111 

3 0.13290 0.037037037 

4 0.05788 0.012345679 

5 0.03646 0.004115226 

6 0.01791 0.001371742 

7 0.01152 0.000457247 

8 0.00515 0.000152416 



TABLE L r^ix as a function of the iteration-order d for LOG 
and TWB. 

and its natural invariant density can be exactly de- 
termined, being expressed in the fashion 

Pznv{x) = , }^ ^ ■ (27) 

TTy/x^l — X) 

LOG is paradigmatic because it is representative 
not only of maps with a quadratic maximum, but 
also emerges when the Lorenz procedure is applied 
to many continuous attractors with basins that may 
be approached with the Lorenz method via a IP- 
map (like the Lorenz, Rossler, Colpits ones among 
others). A non uniform natural invariant density 
is an important feature in this instance [^. The 
ensuing rmia^ -values are also displayed in Table [J 
They have been obtained by recourse to the Trans- 
fer Operator Method, as described in Q. 

• TWB is given by 

{3x„ if < a;„ < 1/3 
3xn - 1 if 1/3 < xn < 2/3 . (28) 
3xn - 2 if 2/3 < Xn < 1 

TWB is representative of the class of piecewise- 
lincar maps as, for example, the Four Way Tailed 
Shift Map, the Skew Tent Map, the Three Way 
Tailed Shift Map, etc. All these maps share a uni- 
form natural invariant density Q . The mixing con- 
stant of the whole family of maps is given by 
<.. = (l/3)' (see table H. 

For the evaluation of the different quantifiers we used 
files with M — 50 ■ 10^ floating point numbers. In the 
Band-Pompe approach we consider D = 6 while for his- 
tograms we have taken Nt,in = 2^^. All Recurrence Plots- 
measures depend on several parameters: 

• the dimension of the embedding space. In this 
paper De = 1. 

• e, a parameter crucial so as to define just when 
recurrences occur. We adopted e — 1/(2^^ — 1) 
corresponding to 16-bits numbers. 



• Imin is the minimum length accepted for diagonal 
lines. Imin = 2 is used in this paper for all diagonal 
measures except for L {Imin = 1 is used for L). 

• N is the number of values used for each realization. 
In this paper values of RR, DET, ENTR and L 
are mean values over 10 surrogate-series with N = 
10000 data each. 

Figures [U [2l and [3] illustrate the behavior of all the 
quantifiers for the iterates of LOG (Figs, (a)) and the 
iterates of TWB (Figs, (b)), respectively. These figures 
show that the following quantifiers are the ones usable 
for measuring r„i^, namely, C^^^^, DET, ENTR, and 
L. On the other hand, the following quantifiers depend 
on the invariant density but they do not depend on rmix ■ 

Intrinsic computation quantifiers display a completely 
different behavior for LOG and for TWB. In LOG both 
quantifiers have no dependence with rmix, but in TWB, 
hf^ decreases as rmix increases while E is an increasing 
function of r^ix- Thus, these parameters do not seem to 
be convenient ones. 

Comparing LOG with TWB by recourse to these pa- 
rameters shows that TWB is slightly better than LOG. 
The problem with TWB and with other piecewise linear 
maps is they are not realistic enough and that their im- 
plementation is more involved that for other maps like 
LOG. 

As an application of the above quantifiers we study two 
usual randomization procedures by means of the repre- 
sentation plane depicted in Fig. 31 employing a quan- 
tifier depending on ^{x) as x-axis 

(^^(hist) selected) 
and a quantifier depending on rmix as y-axis {DET is 
selected). The first procedure is Skipping and the sec- 
ond one is Discretization (see 125||). Skipping has 
been used as a randomization procedure for piecewise 
linear maps. The representation plane evidences the fact 
that this procedure is better than Discretization be- 
cause these maps already have the ideal invariant mea- 
sure (they have ijC****) = 1) and only rmix must be di- 
minished to get the ideal PRNG. Fig. 0] (b) shows that 
Discretization is a better procedure for LOG because 
the ideal point [1,0] is not reached by Skipping. 

V. CONCLUSIONS 

In summary we were able to show here that: 

1. Two classes of quantifiers are required for the eval- 
uation of the quality of a PRNG: a) quantifiers de- 
pending on rmix only (and not on ^{x), like H^'^^'^^\ 
(J {hi St) ^ and RR and b) quantifiers depending on 
IJ,{x) only (and not on rmix), as C^^^\ H^^^\ 
DET, L and ENTR. 

2. Intrinsic computation quantifiers are dependent on 
both /i(x) and rmix and then they are not conve- 
nient for PRNG-analysis with our methodology. 
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3. Representation planes with one quantifier of each 
class as coordinate axis allow for different chaotic 
PRNGs to be compared to each other so as to de- 
termine the best one. Furthermore, these represen- 
tation planes permit one to judiciously select the 
best randomizing procedure. 

Our present results are consistent with those of previ- 
ous works milllil. 
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